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Abstract. We study an extended two-dimensional Shastry-Sutherland model in a magnetic 
field where besides the usual Heisenberg exchanges of the Shastry-Sutherland model 
two additional SU(2) invariant couplings are included. Perturbative continous unitary 
transformations are used to determine the leading order effects of the additional couplings 
on the pure hopping and on the long-range interactions between the triplons which are the most 
relevant terms for small magnetization. We then compare the energy of various magnetization 
plateaux in the classical limit and we discuss the implications for the two-dimensional quantum 
magnet SrCu2(B03)2. 



1. Introduction 

The two-dimensional frustrated quantum magnet SrCu2 (603)2 displays a fascinating sequence 
of magnetization plateaux i2j, which has triggered a lot of activities. Most theoretical studies 
have been devoted to the properties of a 2D spin-1/2 Heisenberg model known as the Shastry- 
Sutherland model[3| in a magnetic field. In the parameter regime relevant for SrCu2(B03)2, 
the ground state of this model is exactly given by the product of dimer singlets [3], and the 
elementary excitations of the model are elementary triplets (triplons |! 5]). The magnetization 
process can then be described in terms of polarized triplons having a magnetic quantum number 
Sz = +1. The polarized triplons are hardcore bosons which interact and move on an effective 
square lattice|6l[7j. The competition between the kinetics and the interaction of the triplons can 
lead to a rich phase diagram since a finite magnetic field creates a finite density of triplons [8]. For 
the case of the Shastry-Sutherland model it is known that the interaction between the triplons 



is the dominant part which is a consequence of the strong frustration. A strong interaction 
leads to the formation of magnetization plateaux which correspond to Mott insulating phases, 

1. e. triplons are frozen in the ground state of the system in a regular pattern. 

All theoretical approaches agree on the presence of magnetization plateaux at 1/3 and 
l/2[6l [71 [U [TOl [TT], in agreement with experiments [H [12]. However, the structure below 1/3 is 
rather controversial. On the experimental side, the original pulsed field data have only detected 
two anomalies interpreted as plateaux at 1/8 and l/4pp, but the presence of additional phase 
transitions and of a broken translational symmetry above the 1 /8 plateau has been established 
by recent torque and NMR measurements up to 31 T[T31 [H]. The possibility of additional 
plateaux has been pointed out by Sebastian et al[T2], who have interpreted their high-field torque 
measurements as evidence for plateaux at 1/q with 2 < q <9 and at 2/9. On the theoretical side, 
the situation is not settled either. The finite clusters available to exact diagonalizations prevent 
reliable predictions for high-commensurability plateaux, and the accuracy of the Chern-Simons 
mean-field approach initiated by Misguich et al. [10] and recently used by Sebastian et al. [12] to 
explain additional plateaux is hard to assess. The essential difficulty lies in the fact that, since 
plateaux come from repulsive interactions between triplons, an accurate determination of the 
low-density, high-commensurability plateaux requires a precise knowledge of the long-range part 
of the interaction. 

The long-range part of the interaction has been determined recently by perturbative 
continuous unitary transformations(PCUTs) |15j . It has been found that only a few two-body 
interactions are decisive at small values of the magnetization and a surprising sequence of 
plateaux at magnetizations 1/9, 2/15, and 1/6 has been deduced by analyzing the effective 
model in the classical limit |15j. This has the advantage that almost all commensurabilities can 
be considered. Note that the effective model obtained by PCUTs is in remarkable accordance 
with a recent contractor renormalization (CORE) calculation[16J. 

The obvious contradictions between current theoretical calculations for the magnetization 
curve in the Shastry-Sutherland model and experimental investigations of the compound 
SrCu2(B03)2 certainly opens the question whether the theoretically considered Shastry- 
Sutherland model is sufficient or if additional magnetic couplings are necessary. A recent ab initio 
calculation^] has suggested possible extensions of the Shastry-Sutherland model which can be 
either SU(2) invariant couplings of longer range (in-plane or intra-plane) or Dzyaloshinskii- 
Moriya (DM) interactions. In this work we will study the leading order effects of the two 
dominant SU(2) invariant in-plane couplings. 

2. Model 

The total Hamiltonian studied can be written as 

H = Hss + H, + Hi-BY^St , (1) 

i 

where ilfss denotes the usual Shastry-Sutherland model with Heisenberg-type couplings J and 
J' (see Fig. [1^). The first extra term represents a Heisenberg-type coupling oc Jc which gives 
rise to dimerized chains in Jc along both diagonals x it y of the Shastry-Sutherland lattice (see 
Fig. [T}3). Similarly, the second extra term H\ is also a Heisenberg-type coupling oc J\ which 
connects dimers along both diagonals x ±y, but this time two-leg ladders with a rung coupling 
J and a leg coupling Ji are formed (see Fig. [TJ;). Lastly, the term proportional to B represents 
an external magnetic field in z-direction. 

We are actually not interested in the full phase diagram of the model. A clear hierarchy 
is expected for a theoretical description of SrCu2(B03)2: J > J' ^ Jc, Jy- Recent ab-initio 
calculations estimated that Jc ~ 0.025 and Ji «0.01 J[17j. One may wonder why such small 
couplings can have a sizable effect on the magnetization process. To this end it is important 



to realize that at low magnetization the relevant two-body interactions and the pure hopping 
start only in sixth order in J'/J[T5] which is also small quantity. One finds (J'/J)^ = 0.015625 
for J' / J = 1/2 which is indeed of the same order of magnitude. It is therefore mandatory to 
perform a high-order expansion in J'/ J, but it is enough to consider only effects which are linear 
in Jc and J\ because higher orders are negligible. The first part has been done in Ref. [15J. In 
this work we determine the corrections of the form J^PdJ' / J) and J\P\{J' / J) where Pq and Pi 
are polynoms in J' / J which we have determined up to order 5. 




Figure 1. Figure illustrates the different couplings studied in this work: (a) usual couplings 
J and J' of the Shastry-Sutherland model, (b) coupling creating dimerized chains along 
the diagonals, and (c) Ji creating two-leg ladders along the diagonals. 



3. Method 

We have tackled Hamiltonian ([1]) by PCUTs as we did recently for the pure Shastry-Sutherland 
model in a magnetic field[T5]. We refer for a detailed introduction of the method to the existing 
literature [181 [T9| I20j. Here we only give the general idea and we stress certain technical aspects 
which arise due to the new couplings. 

The main idea of the PCUT approach is to transform the initial Hamiltonian ([T]) which 
changes the number of triplets to an effective Hamiltonian which conserves the number of the 
true quasi-particles which are triplons in our case. In PCUTs this is done efficiently to high order 
in perturbation. It is important to realize that the magnetic field term is unchanged during the 
unitary transformation because the total S^ot is a conserved quantity. The relevant processes 
for the physics in a finite magnetic field have maximum total spin and total Sz- Other spin 
channels are relevant for spectroscopic observables which have been studied earlier for the pure 
Shastry-Sutherland model [211 [22]. The effective Hamiltonian for the physics in a finite magnetic 
field therefore contains only one triplon flavor with 5^ = 1. This remaining degree of freedom is 
a hardcore boson [BJ [151 120] . 

There are several reasons why the treatment of Hamiltonian ([T|) is more complicated than 
the pure Shastry-Sutherland model. First of all, there are three expansion parameters {J' / J , 
Jc/J-, and J\/ J) compared to one {J' /J). Second and more important, the product state of 
singlets is not the exact groundstate anymore once Jc and/or Ji are finite. This results in much 
more intermediate states during the numerical calculation and leads consequently to a reduced 
maximum order which can be reached by PCUTs. Finally, the number of matrix elements 
of Hamiltonian ([T|) are larger than in the pure Shastry-Sutherland model, e.g. there are only 
processes changing the number of triplets by one in the Shastry-Sutherland model while He and 
Hi contain both matrix elements changing the number of triplets by two. This results also in 
an increased effort for the PCUT calculation. 



4. Effective couplings 

In this section we will present the results we have obtained for the most important processes 
at low magnetizations. These are on the one hand the pure kinetic hopping which is important 
for possible superfluid and supersolid phases and on the other hand the two-body interactions 
which determine the structure of the magnetization plateaux. 

4.I. Kinetics 

For the pure Shastry-Sutherland model, a direct hopping of the triplon is heavily suppressed due 
to the strong frustration. The kinetic energy is dominated by correlated hopping processes [B] 
which favor so-called supersolid phases[23J. This results in an almost flat one-triplon dispersion. 
The only non-vanishing hopping in the effective model up to order 15 is the hopping over 
the diagonals in the Shastry-Sutherland lattice (see Fig. As mentioned above, it only starts 
in sixth order with a very small coefficient (l/96)( J'/J)^. Taking the full expression, one gets 
tf/J = 0.00045 for J' /J = 0.5 which is a very small energy scale. 




Figure 2. Shastry-Sutherland lattice and definition of the 2-body interactions. Vn is the 
coefficient of the 2-body interactions between the thick dimer and the dimer labeled Vn- 
Additionally, the definition of the hopping t2 is illustrated in blue on the right. 

It is then obvious that the couplings Jc and Ji will dominate the pure kinetics of the triplons 
even if they are only a few percent of J. The reason is that they are contributing to t2 linearly in 
Jc and J\. As mentioned above, the term He represents dimerized chains along both diagonals of 
the Shastry-Sutherland lattice. The corresponding matrix element for a hopping of the triplon 
over the diagonal is in first order t2 = — Jc/4. In contrast, the term H\ forms two-leg ladders 
over the diagonal of the Shastry-Sutherland lattice and one finds in first order = Ji/2. Note 
that the different signs of the two contributions arise due to the frustration. Indeed, the coupling 
Jc is a frustrating diagonal coupling in the two- leg ladders formed by the couplings Ji. 

In total, we obtain the following hopping over the diagonal 

Assuming J^, » Jj » O.OIJ, we obtain J = 0.00295 for J'/ J = 1/2. The pure hopping has 
therefore increased by a factor 10 due to the extra couplings! Note that there will be no sizable 
effects on other hopping integrals. 



4-2. Two-body interactions 

In Ref. [15], we have studied the effective model obtained by PCUTs by analyzing the classical 
limit, i.e. we replaced the hardcore bosons by effective spin 1/2 using the Matsubara-Matsuda 
representation [24] of hard-core bosons, and we treated the spins as classical vectors of length 1/2. 
We found that the magnetization curve is dominated by solid phases. Each solid corresponds 
to a phase with broken translational symmetry where triplons are frozen in a regular fashion in 
the ground state. These gapped phases produce plateaux in the magnetization curve. 
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Table 1. This table illustrates the effect of Jc and J\ on all the considered two-triplon 
interactions. For a specific interaction Vn (column 1), the leading order concerning the three 
expansion parameters J', Jc, and Ji are given in columns 2-4. Column 5 shows the amplitude 
of the two-triplon interaction in the Shastry-Sutherland model for J'/ J = 1/2 using the 
extrapolated 15*^^ order series. Column 6 (column 7) gives the amplitude of the bare series 
linear in Jc (Ji) for J'/J = 1/2 and Jc/J = 1/100 (Ji/J = 1/100). 

Surprisingly, the plateaux stabilized at low magnetizations are 1/9, 2/15, and l/6|15j which 
is in contradiction to the expected 1/8 plateaux in the experimental compound SrCu2(B03)2. 
The classical energy of the different plateaux mainly depends on the two-triplon density-density 
interactions as long as we consider low magnetization. The different two-triplon interactions are 
illustrated in Fig. [2j One central result of Ref. [15j is that the classical energy of all stabilized 
structures at low magnetization (except the 1/9 plateau) depends dominantly on the three 
interactions Kj', V5, and V^ which all appear in sixth order (J'/J)^. This statement is also true 
for the 1/8 structures discussed in the literature |llj. We will therefore discuss in the following 
the influence of Jc and Ji on all the considered two-triplon interactions but with a special focus 
on the interactions V^, V5, and Vj. 

In Tab. [1] we give an overview of the effects of Jc and Ji on the two-triplon interactions. 
As expected, the largest interactions are not much affected by the additional couplings. These 
terms appear in low order in J'/J and the change is only a small perturbation. But this is 
different for the interactions starting in higher order in J'/J which are the relevant ones at 
low magnetization. On the right side of Tab. [1] we give a quantitative comparison between the 
interaction strength in the pure Shastry-Sutherland model at J'/J = 1/2 (using the extrapolated 
25th Qpjjgj. series) and the additional contributions linear in Jc and J\ setting each coupling to 
O.OIJ. Note that this is indeed the order of magnitude deduced by ab-initio calculations for 
these couplings (in fact one obtains slightly larger couplings of the order 0.025 J for Jn)[T7|. 

Among the relevant interactions for the physics at low magnetization, the coupling affected 
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Table 2. The left table shows the classical energy of the relevant plateaux at low magnetization 
in terms of the twp-triplon interactions. The right figure illustrates the different plateau 
structures and their unit cells. Black dimers denote triplets and empty dimers represent singlets. 



most by Jc and Ji is the two-triplon interaction Vj which gets already a contribution of the order 
Jc{J' / J)"^ and J\{J' /J)"^. This is related to the fact that the interaction V3 starts already in 
order (J'/J)^. Indeed, the distance between the two dimers for the Vy interaction is the distance 
for the V3 interaction plus one diagonal link. The diagonal link can be done linear in J^, 1 such 
that one finally has Jc^\{J'/J)^. 

Globally, one also recognizes that the effects of Jc and Ji on the two-triplon interactions are 
often competing, i.e. although each change can be sizable, the combined effect of both couplings 
almost compensate assuming Jc ~ J]. 

5. Magnetization plateaux 

In this section we describe how the additional couplings Jc and Ji affect the classical energy 
of the relevant plateaux and its structures at low magnetization. Here we concentrate on the 
three structures at magnetizations m G {1/9,2/15, 1/6} we found to be stabilized for the pure 
Shastry-Sutherland model[T5] plus two structures for 1/8 plateaux discussed in the literature for 
SrCu2(B03)2[9][l] The explicit expressions for the classical energy of these plateaux in terms of 
the two-triplon interactions Vn are given in Tab. El Additionally, the corresponding structures 
of the plateaux are illustrated. 

In a next step we use the values for the two-triplon interactions as discussed in the last section 
in order to compare the energy of the different plateaux. This is done in Fig. [3] for different 
representative values of J', Jc and Ji. For the pure Shastry-Sutherland model (see Fig. [3ti), the 
sequence 1/9, 2/15, and 1/6 is stable. 

Adding only the coupling Jc, results in an increase of all the couplings V^, V5, and Vj. 
Correspondingly, all classical energies are larger and the phase transitions occur at larger values 
of the chemical potential (see Fig. EJ)). The coupling affected most is the two-triplon interaction 

^ Note that we do not discuss in this work the possibihty that the new couphngs Jc and Ji induce novel solid 
structures or stabilize suprafluid or supersolid phases. We concentrate here on the plateaux stabilized for the 
pure Shastry-Sutherland model plus additional structures at 1/8 magnetization. 




Figure 3. Classical energy of the different low-magnetization plateaux for J' / J = 1/2 and (a) 
J^ = Jy = 0, (b) Jc = 0.05 and Ji = 0, (c) Jc = and Ji = 0.05, and (d) = Jj = 0.05. 



Vi which affects mostly the 1/6 plateau, then the 2/15 plateau and also the 1/8 rhomboid 
structure. One therefore obtains a larger 1/9 plateau and also the energy of the 1/8 square 
structure has only a slightly larger classical energy close to the transition from 1/9 to 2/15. The 
Jc coupling therefore tends to stabilize a 1/8 plateau. For Jc > 0.05, we indeed find a stable 1/8 
plateaux with a square unit cell. Note that the 2/15 plateau has an increased extension because 
the energy of 1/6 structure has changed even more. 

We next consider the case where only the coupling J\ is turned on. This is illustrated in 
Fig. [3b. Here the interactions and V5 are almost unchanged. The interaction Vj is again 
strongly altered but this time it gets a negative contribution. So the effects are opposite to the 
coupling Jc. The 1/6 structure benefits most and its extension strongly increases. We would 
like to remark that the case of a finite J\ with vanishing Jc seems to be unrealistic since it is an 
exchange over a longer distance than Jc. 

Finally, we consider the case where both couplings are present having the same weight. A 
corresponding plot is shown in Fig. [3li. It should be clear from the discussion above that the 
effects on the interaction V7 are almost compensating. One mainly has an increase of and V5 
due to the coupling Jc. Since all plateaux structures except 1/9 depend on these two interactions, 
the net effect is rather small. Only the 1/9 plateau has an increased extension. 

6. Conclusions 

We have studied the magnetization process of an extended Shastry-Sutherland model. Two 
SU(2) invariant in-plane couplings have been added to the normal Shastry-Sutherland model. 
Perturbative continuous unitary transformations are used to derive the leading order effects of 
these couplings on the relevant hopping and interaction amplitudes. 

The most important finding of this work is that already very small couplings can affect the 



physics at low magnetization significantly. The physical reason is that the relevant processes in 
the effective model for the pure Shastry-Sutherland model in this regime start only in 6*^ order 
perturbation theory which is also a small energy scale. We found that the additional coupling 
which is expected to be the largest in-plane correction to the Shastry-Sutherland model indeed 
stabilizes a 1/8 plateau with a squared structure but the coupling strengths needed for a sizable 
1/8 plateaux are a factor 3 to 4 larger than the value deduced from first principle studies [T7]. 
In contrast, the sum of both couplings assuming Jc ~ J\ has only a small net effect which is a 
consequence of their compensating character. 

Our study clearly shows that it will be also interesting to investigate the effects of 
other additional couplings to the Shastry-Sutherland model like intra-plane couplings or DM 
interactions on the magnetization process which are expected to be of similar magnitude than Jc 
and J\ in order to gain a better understanding of the rich magnetization curve observed 

for the frustrated quantum magnet SrCu2(B03)2. Additionally, it will be also important to 
study the effect of the enhanced kinetic energy induced by the additional couplings in order to 
investigate whether superfluid or supersolid phases are formed. 
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